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Abstract 

In Boundary Element Method, Green's function with no boundary 

conditions is used for solving Laplace's equation with Dirichlet boundary 
condition. To determine the gradient of solution on the boundary, we 
need to solve the boundary integral equation numerically in most practi- 
cal cases. 

Here we discuss the alternative method to avert solving that boundary 
integral equation. It is based on the solution for Poisson's equation which 
has the singularity on the boundary specified by the boundary function 
and shown it is applicable to arbitrary shaped domain and reduces calcu- 
lation cost considerably. 

keywords Boundary Element Method, Dirichlet problem, Poisson's equa- 
tion. Green's function, boundary integral equation 



Introduction 

In Boundary Element Method (BEM), Laplace's equation is solved by us- 
ing the following Green's theorem; 

{(t)V^i:-ijV^(t))dV = [ {<j)Vi)-i)V(j))da (1) 

D JdD 

Let V^(/)(a;) = and ^ the Green's function G for the Laplacian with no bound- 
ary conditions. Then, 



4){x)= [G{x,x')V<t){x')- (t){x')VG{x,x')]da{x') (2) 

JdD 

where x is in the domain D. When x is on the boundary dD and dD is suffi- 
ciently smooth around x, (j){x) satisfies the following equation; 

\(j>{x)= [ [G{x,x')V<p{x') - (l){x')VG{x,x')]da{x') (3) 

2 JdD 



Under Dirichlct boundary condition, it is necessary to solve the Fredholm inte- 
gral equation of the second kind shown in (3) for determining the gradient of 
solution on the boundaryj^l 

The boundary integral equation (3) is not solved analytically unless the 
boundary has some special forms [''1 and usually solved numerically for practical 
purposes. It requires to deal with the large and fully populated matrix repre- 
senting that equation for estimating the solution with sufficient precision, fur- 
ther the singularity caused from Green's function must be treated carefully J^l 
Therefore generally the numerical solution of boundary integral equation re- 
quires great computational cost, and it is expected to improve the performance 
of BEM considerably by replacing this process with any alternative methods 
which don't require estimating the gradient of solution on the boundary. In the 
following section we discuss how to construct the method for satisfying such 
request. 



Second partial derivatives of boundary function for solving Dirichlet 
problem 

One of the methods to solve Dirichlet problem without estimating the gradi- 
ent of solution on the boundary uses Poisson kernel. For example, the Dirichlet 
problem for the unit disk in the complex plane is solved by using Poisson kernel 
as, 

where f{z) is the boundary function. Using this formula, the Dirichlet problems 
for various regions in the complex plane are solved by mapping unit disk to the 
region conformally^l or applying the balayage method to the sequence of disks 
approximating the region. 1^1 However practically constructing such conformal 
map or applying the balayage method requires considerable effort when that 
region has complex shape. Here we discuss not only Poisson kernel, but also 
alternative boundary function in order to solve Dirichlet problem. 

The above alternative boundary function is actually considered on the Pois- 
son's equation V^^ = g where the solution <j) coincides with that of Laplace's 
equation V^V = satisfying Dirichlet boundary condition: tp = f and cj) = tp in 
the domain surrounded by arbitrary shaped but sufficiently smooth boundary. 
Namely we are to look for the function g which makes (j) = ip m that domain. 

Firstly we discuss two-dimensional case in real plane R^. Let f{x,y) the 
Dirichlet boundary condition of Laplace's equation V'^ip = 0, i.e. il>{x,y) = 
f{x,y) on the boundary dD where dD is sufficiently smooth. We suppose D is 
closed and the boundary function is sufficiently smooth. For each point (x, y) 
on the boundary dD we take the local coordinate system {x',y') where x' is 
tangential on the boundary point and y' is orthogonal to it. The origin of that 
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coordinate system is the boundary point {x,y). Also we suppose the Jacobian 
§^ is equal with 1 (Fig.l). 




Fig.l Local coordinate 
f^^\x,y) is defined on the point r = (x,y) e dD as 

dx 



(x',y')=(o,o) {{x' ,y') -.local coordinate) (4) 



We suppose f^^\x,y) is also sufficiently smooth. Using f^^\x,y), f^'^\x,y) is 
defined as, 



}^^\x,y) = PHx.y)5{y') 



mi 



(5) 



where 5{y') is one-dimensional delta function. Then the following proposition 
is satisfied; 

Proposition 1. Let xl) the solution of Laplace's equation: 

= (6) 

under the Dirichlet boundary condition ip{x,y) = f{x,y) on dD, then there 
exists (j) such that, 

ip = (p + c (c = constant) : (a;, y) G D 

where <p satisfies the following Poisson's equation; 

VV(a;,y) = f^''\x,y) : {x,y) € OD or 



= 



■.{x,y)^dD 



(7) 
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Proof. From the Poisson's equation (7), the Laplacian of (j) at the point 
{x,y) on the boundary is f^^\x,y). Since f'^^\x,y) includes delta function, the 
Laplacian diverges on the boundary. However it does not mean that the second 
partial differential of d along every direction always diverges on the bound- 
ary. To see it, firstly we define the function g{x,y) as satisfying the following 
condition in the neighborhood of (x, y). 

y^g{x, y) = and 

dM^',y') , _ dy{x',y') ^ 

l(x',i/')=(0,0) - \(x',y') = {0,0) 

where {x',y') is the local coordinate around the point {x,y), then a solution (p 
of Poisson's equation (7) is expressed in the neighborhood of {x, y) as, 

(l>{x, y) = g{x', y') -.if y' <Qor 

d^f{x',y'), 



dx'"^ 



\(x' ,y')=(ofi) ■ y' + 9{x' ,y') ■.ify'>0 (9) 



where the origin of the local coordinate is taken on dD. From the above ex- 
pression, wo know that the second partial differential of (f) along x' at the point 
{x, y) on dD is finite, i.e. 

l(x',y') = (0,0) - \(x',y') = (Qfi) U^j 

Therefore cj) is expressed in the neighborhood of (x, y) fixing y' = as, 

(j){x' , 0) = f{x', 0) -|- c\x' + C2 (ci, C2 : constants) (11) 

Owing to the periodicity of ((>{x,y) and f{x,y) on dD, the constant ci in the 
above equation vanishes. On the other hand the solution 'ip of Laplace's equa- 
tion: 

vV = o 

under the Dirichlet boundary condition: t/'(a;, y) = f{x, y) + C2 is expressed with 
ip in (6) as, 

V' = + C2 (12) 

Owing to the uniqueness of the solution of Dirichlet problem, ^ coincides with 
(j) in the domain D, i.e. 

(f> = i} + C2 (13) 

in D. II 
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Since /*^^^ (r) is the product of a smooth function and one-dimcnsional deha 
function as defined at (5), the solution </> of the Poisson's equation (7) expressed 
with Green's function as, 

0(r) = / /(2)(r')G(r,r')da' (14) 

JD 

is actually the integral along the boundary dD. Therefore it is converted to 
one-dimensional integral using f^'^\r') in (4) as, 

<P{r) = C [ P\v')G{r,r')dl' (15) 

JdD 

where C is a constant (unit: length). If </> under the condition: (j){r) — >■ when 
r — )• oo is the solution, it must satisfy the following integral equation; 

\(x'.y) = {o,o) 

= ^"^'^ W = L /^-HrO ^'^"^:l"''\ ^'.')Ko,o) dl' (16) 

where r is taken on dD. Namely — ^ is the eigenvalue of the above integral 
equation. As discussed later, C is estimated numerically without solving that 
eigenvalue problem in the application to BEM. 

The second partial differential of the boundary function f{x,y) shown in 
(4) is actually carried out in the tangent space Tp{dD) of each point on the 
boundary dD. Therefore when the boundary is not smooth enough to define 
the tangent space, the above proposition is not satisfied. 

The above proposition is extended to higher dimensional one with the man- 
ner similar to two-dimensional case. In the n-dimensional case we suppose the 

boundary dD is compact. The local coordinate {x[, ...,x'^_i,Xj^) at the point 
p on the boundary is set so as to make orthogonal to the tangent space 
Tp{dD) and ^l^^' ''^'') = 1- Also its origin is p. /'•^•'(r) is defined on the point 
r = (xi, ...,x„) e dD as, 

J nr) - 2^ |(xi,...,<)=(o,...,o) (17) 

and /'•^•'(r) is defined as, 

f'\r) = p\r)S{x'J (18) 

where S{x'^) is one-dimensional delta function. Then the following proposition 
is satisfied; 

Proposition 2. Let tp the solution of Laplace's equation: 

VV = (19) 
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under the Dirichlet boundary condition il){r) = /(r) on dD, then there exists ^ 
such that, 

ip = (f) + c (c = constant) : {x, y) & D (20) 

where (p satisfies the following Poisson's equation; 

VV(r) = f^^\r) :r€dDor 

= ir^dD (21) 



Proof. The function g{r) is defined as satisfying the following condition 
in the neighborhood of the point r = {xi,...,Xn) in R"; 

V^5(r) = and 

E""^ d'^g{x[,...,x'J 
\(x[,...,x'J={0,...,0) 

j=i j 

E^^/(^'li---'^n)| /oo\ 
IK,...X) = (0,-,0) {^^) 

j=l j 

Then a solution </) of Poisson's equation (21) is expressed in the neighborhood 
of (a;i, ...,Xn) as, 

0(r) = g{x[,...,x'J :ifx'^<Oor 

_ ''y^ d^f{x[,...,x'J 

~ 1^ Qx^ IK,...,x;)=(o,...,o) • 

+ :i/a;;>0 (23) 

where the origin of the local coordinate is taken on dD. It satisfies the following 
equation; 

d'^<j){x[,...,x'J ^ _'i^ d^f{x[,...,x'J ^ 
^ dx^ IK,...,x;.)=(o,...,o) - 2^ l(x;,...,x;)=(o,...,o) [^V 

at the points on the boundary dD. Set cr(r) = (f){v) — /(r), then the equation 
(24) is rewritten as a locally expressed Laplace's equation, i.e. 

VV(r) 

dx'l dx'^_, 
= (25) 

in the tangent space of each point r on the boundary dD. 
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As long as (T(r) is sufficiently smooth, (T(r) docs not take maximum (or 
minimum) value in dD, otherwise at the maximum (or minimum) point the 
Laplacian V'^a{r) would not be equal with zero. Since dD is compact, 0-(r) is 
a constant on dD, i.e. 

0(r) — /(r) = c (c : constant) (26) 
On the other hand the solution ^ of Laplace's equation: 

V^ = 

under the Dirichlet boundary condition ip{r) = /(r) + c is expressed with ip in 
(20) as, 

= V + c (27) 

Owing to the uniqueness of the solution of Dirichlet problem, 'ip coincides with 
(j) in the domain D, i.e. 

</> = V + c (28) 

in D. II 

If the boundary dD is not compact, the solution of Laplace's equation (25) 
may take the maximum (or minimum) value on the closure of dD. Therefore it 
may not be constant and the proposition 2 is not satisfied in general. 

The solution (p of the Poisson's equation (21) expressed with n-dimensional 
integral as, 

0(r)= / /(2)(r')G(r,r')dV^' (29) 

JD 

is converted to (n-l)-dimensional integral as well as the two-dimensional case 
discussed in the proposition 1 as, 

<P{r) = C [ P\v')G(v,v')da' (30) 

JdD 

where C is a constant (unit: length). Actually C is the eigenvalue of the fol- 
lowing integral equation; 




,...,0=(o,...,o) 



= P\r) = cY^[ /(^HrO^^^^%^lK,...,<)=(o,....o)M31) 

j=l J9D OX- 

where r is taken on dD. 
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Application to BEM 

Prom Proposition 2, tlie solution of Laplace's equation under Dirichlet bound- 
ary condition is expressed in general with the boundary integral (30) as, 

V(r) =C\ [ P^G{r, r') da' + C2 (32) 

JdD 

Since both /'^^^(r) and G(r,r') arc explicitly calculable, the V'(r) expressed with 
(32) is estimated without any unknown factors except constants Ci and C2. 
Once we have got Ci and C2, we don't have to estimate the gradient of solution 
on the boundary, and it means the process for solving the boundary integral 
equation (3) has been eliminated. 

In order to apply the equation (32) to BEM, the boundary and boundary 
function must be finite and sufficiently smooth. The techniques for constructing 
mesh on the boundary used for ordinary BEM is effective as well. As the simplest 
way of numerical evaluation, the equation (32) is discretized as, 

JV 

Vi(r) = Ci ^p)(r,)G(r,r,)Aa(r,) + C2 (33) 

where /^^^ is defined by discretizing (17) as, 

f(2). ^ AV(a;;,...,<) Ay{x[,...,x'J A^f{x[,...,x'J 
^ Axf + Ax'i Ax'^_, 

Ci and C2 are obtained by solving the following linear system; 

N 

i>{r,) = Ci5]/(2)(r,)G(r„,r,-)Aa(r,-) + C2 (35) 

TV 

= C,J2PHrj)G{n,rj)Aa{rj) + C2 (36) 

where Va and ri, are the mesh points chosen arbitrarily on the boundary. V'(''a) 
and ip{''^b) are given previously as the boundary condition. 

Since the calculation loads for obtaining each f^^\rj), G(r, r^), and Acr(rj) 
don't depend on the number of mesh points A'^, the order of calculation load for 
estimating ^{r) inside of the boundary is 0{N). 

In the Pig. 2, two-dimensional Laplace's equation solved by the above method 
is illustrated. The Laplace's equation dealt with here has the following boundary 
conditions; 

dD = {{x,y)\x'+y^ = l} 
f{cos0,sine) = sin2e {0 < 9 < 2Tr) 
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Fig. 2 Numerical evaluation of Poisson's equation having singularity along the 

boundary 



Actually it is the numerical solution of Poisson's equation (7) and defined not 
only in the domain D for the Laplace's equation but also outside of it. As 
seen in the figure, the second partial differentials of the solution diverge on the 
boundary except the direction tangential to the boundary as mentioned in the 
proof of Proposition 1. 



Summary 

The solution of Laplace's equation under Dirichlet boundary condition co- 
incides with the solution of a Poisson's equation which has a specific singularity 
along the boundary. The domain in which the solution of that Poisson's equa- 
tion is applicable is arbitrary as far as the boundary of the domain is compact 
and the boundary function is finite and both of them are sufficiently smooth. 

Comparing the solution derived from the above Poisson's equation with that 
from ordinary Boundary Element Method, we find the former does not need to 
solve the boundary integral equation which is used for estimating the gradient 
of the solution on the boundary and requires considerable calculation cost in 
the latter. Therefore it is expected to reduce the calculation cost by implement- 
ing the method discussed here to Boundary Element Method especially when 
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dealing with the Laplace's equation in the domain having complex shape. 
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